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ABSTRACT 

We have investigated the accuracy of the Is-vacancy fluorescence data base 
of Kaastra & Mewe (1993) resulting from the initial atomic physics calculations 
and the subsequent scaling along isoelectronic sequences. In particular, we have 
focused on the relatively simple Be-like and F-like Is-vacancy sequences. We 
find that the earlier atomic physics calculations for the oscillator strengths and 
autoionization rates of singly- charged B II and Ne II are in sufficient agreement 
with our present calculations. However, the substantial charge dependence of 
these quantities along each isoelectronic sequence, the incorrect configuration 
averaging used for B II, and the neglect of spin-orbit effects (which become 
important at high-Z) all cast doubt on the reliability of the Kaastra & Mewe 
(1993) data for application to plasma modeling. 

Subject headings: atomic data - atomic processes - line: formation - X-rays: 
general 
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1. Introduction 

In collisionally ionized or X-ray photoionized plasmas, high-energy electrons or photons 
lead to the production of Is-vacancy ionic states which then decay via sequential emission 
of single or multiple electrons and/or photons. The exact strengths of these competing 
processes determine fundamentally important quantities of the plasma such as the ionization 
balance and the observed spectra of emitted and/or absorbed photons. Hence, interpreting 
the properties of these plasmas requires accurate atomic physics calculations for the various 
autoionization and radiative rates. Here we are interested in assessing the accuracy of the 
available data base that provides such computed (or inferred) Auger rates and fluorescence 
yields to the astrophysics community. The accuracy of these atomic data are crucial to the 
interpretation of the spectra of photoionized plasmas such as are found in X-ray binaries and 
active galactic nuclei. These data are also important for supernova remnants (SNRs) under 
conditions of non-equilibrium ionization (NEI). 

Two of the more widely used spectral codes for modeling photoionized plasmas are 
CLOUDY (Ferland et al. 1998) and XSTAR (Kallman & Bautista 2001). A commonly used 
code for modeling NEI in SNRs is that of Borkowski, Lyerly, & Reynolds (2001). These all 
in turn rely on the table of electron and photon emission probabilities compiled by Kaastra 
& Mewe (1993). This comprehensive data base considers the sequential multiple electron 
and/or photon ejections for all stages of all Is- vacancy ions in the periodic table up through 
zinc. In order to produce such a massive array of numbers, however, certain approximations, 
questionable from a purely theoretical atomic physics standpoint, were invoked. First, the 
only rigorously computed atomic rates were taken from the early works of McGuire (1969, 
1970, 1971, 1972) for singly- charged ions, which furthermore neglected configuration inter- 
action (CI) and spin-orbit effects. Due to the limited computational resources available at 
the time, and the approximations thus needed to perform such calculations, even these can- 
not be considered as rehable as those that can be carried out with today's state-of-the-art 
capabilities. Second, these singly-ionized results were then scaled along entire isoelectronic 
sequences, assuming constant autoionization rates and oscillator strengths; this approxima- 
tion is least valid for near-neutrals. 

A third approximation used by Kaastra & Mewe (1993) is that the electron and photon 
emission yields were computed using radiative and autoionization rates that were configura- 
tion averaged over possible terms, and the fiuorescence yield was then given as a ratio of the 
averaged radiative rate to the sum of the averaged radiative and averaged Auger rate. For 
modeling purposes, however, this is incorrect; the actual required value is the average of the 
term-specific yields - an average of ratios rather than a ratio of averages. In other words, the 
relative probability of producing each specific inner-shell-vacancy term, and its subsequent 
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term-specific decay, needs to be considered, and this was not done correctly for the data 
compiled by Kaastra & Mewe (1993) (see also Chen et al. (1985) for a further discussion of 
this importance). 

In this paper, we investigate the validity of the above three approximations in order to 
assess the accuracy of the resultant data base of Kaastra & Mewe (1993). To this end, we 
first study the simplest Is- vacancy system that can radiate via a 2p ^ Is dipole-allowed 
transition. This is the removal of a Is electron from the groimd-state B-like sequence, or 
rather the ls2s^2p Be-like inner-shell excited sequence, which is investigated in the next 
section, and which is further simplified by the fact that only one electron, or one photon, 
can be emitted. We follow in Section 3 with a study of the simplest closed- (outer) shell case 
of F-like ions, corresponding to Is vacancies from the Ne-like sequence. A summary of our 
findings and concluding remarks are then given in Section 4. 



2. Case Study of the Be-Like Fluorescence Yields 

Inner-shell Is vacancy of a Be-like ion, whether by photoionization or electron-impact 
ionization of B-like ions (or by photoexcitation or electron-impact excitation of Be-like ions), 
results in either the ls2s^2p(^P) state or the ls2s^2p(^P) state. Prom an independent 
particle perspective, in LS coupling, the following competing decay processes can then occur: 

ls^2s\^S)+uj , (1) 
ls'2sCS)epC''P) , (2) 
ls^2pCP)esC'^P) , (3) 

that is, the Is-vacancy state can either fluoresce, if it is in the state, with a radiative 
rate Ar, or autoionize, from either state, with a total state-dependent rate Aa — Aai + Aa2, 
yielding free electrons denoted by el. (If left in the ^P state the ion does not fluoresce - we 
consider CI and spin-orbit effects in the next section). The radiative rate Ar in atomic units 
(1 a.u. = 4.1341 X 10^^ s~^) is related to the dimensionless emission oscillator strength / by 

Ar = 2u;Wf , (4) 

where u is the emitted photon energy in a.u. (1 a.u. of energy = 27.211 eV) and a ~ 1/137 
is the fine structure constant. Here we define the emission oscillator strength as the absolute 
value of the oscillator strength from the upper ls2s^2p(^P) term j to the lower ls'^2s'^(^S) 
term i: 



1s2s^2pCP) ^ 
ls2s'^2pC^^P) ±1 
ls2s^2p(^'2p) 



/ ^ \fji 



(5) 
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and can thus be related to the absorption oscillator strength fij from the lower term i to the 
upper term j via 

9ifij = 9j \fji\ , (6) 

where Qi = 1 and gj = 3 are the statistical weights of the initial and final Be-like terms, 
respectively. 

Oscillator strengths are more convenient quantities to use along isoelectronic sequences 
because they exhibit certain bounds. Since the absorption oscillator strength is bounded by 
< fij < Ni, where A^j = 2 is the number of Is electrons, the emission oscillator strength is 
bounded by < / < {gi/ gj)Ni = 2/3 (for the present cases), and is a well-behaved function 
of the nuclear charge Z. In fact, if the hydrogenic approximation is valid, i.e., if the nuclear 
potential dominates over the interelectronic repulsive potential, then the emission oscillator 
strength is independent of Z, and the same is true for the autoionization rate Aa- Such an 
approximation is valid for highly-charged ions but not for lower-charged species. 

The fluorescence yield ^, from a given inner-shell vacancy state, is a measure of the 
relative probabilities of the radiative and autoionization decay pathways and is defined as 



Aa 
f 



Thus it only depends on the squared transition energy u"^ and the ratio of the autoionization 
rate to the emission oscillator strength Aa/ f . In the hydrogenic approximation, these scale 
respectively with nuclear charge as and (i.e., independent of g), where g = Z — 3 is the 
asymptotic ionic charge seen by the outer-most electron of the Be-like ion (Cowan 1981). 
With these scaling properties, the expected behaviors at \ow-Z and high-Z are ^ ~ and 
^ 1 (provided / 7^ 0) , respectively. 



2.1. Initial Populations, Configuration Interaction, and Spin-Orbit Effects 

As pointed in Section 2, both the and terms can be populated after Is pho- 
toionization or electron- impact ionization. Following Cowan (1981), and using the sudden 
approximation, we have determined that the probability of populating each term can be 
deduced by considering the squared recoupling coefficient 



\{[{lsls) e^)] 2p (^P) I [(ls2p) ({^^+^>P)] Is {^P))\' = 



(-1)^+^(25+1)^!'^ f f 
25 + 1 



(8) 



- 5 - 

where 5 = for the state and 5 = 1 for the state. This means that the states 
are populated according to their statistical weights, and the ^P state is populated with a 
probability of 1/4. (In general, there also should be a recoupling coefficient involving the 
orbital angular momenta of the three electrons in Eq. 8; however, the / = values for two 
of the electrons' orbital momenta reduces the coefficient to unity for the present case.) We 
have also verified this computationally by performing R-matrix photoionization calculations 
using the Wigner-Eisenbud R-matrix method (Burke & Berrington 1993; Berrington et al. 
1995). Using both approaches we find that in intermediate couphng, the states are also 
populated according to their statistical weights (a similar expression to Eq. 8 involving the 
total angular momentum values j for each electron can be obtained). 

Considering the relative populations of the ls2s^2p Is- vacancy states, the desired quan- 
tity for plasma modeling purposes is the configuration-average fluorescence yield. If CI and 
spin-orbit effects are neglected, this can be defined as an average over LS single-configuration 
(SC) terms as 

_ 1 

z— >oo ^ ' 

where fluorescence from the ^P state is zero so that the asymptotic behavior at large Z is 
1/4. 

CI and spin-orbit effects modify this behavior, however. The largest CI effect is the 
intrashell mixing Cils2s^2p + C2ls2j9^, where the mixing fraction |c2/cip is essentially term 
independent and Z independent for nonrelativistic calculations - it varies between 0.067 
for B II and 0.053 for Zn XXVII. This mixing affects the computed emission oscillator 
strength / and autoionization rate Aa at the near neutral end of the sequence, but changes 
the high-Z ffuorescence yield by less than 10%. The more important CI effect is that the 
admixture of the ls2p^ configuration in the ^P term allows it to radiate to the ls^2p^(^P) 
state. This C2ls2p^(^P) ls^2p^(^P) radiative rate is about a factor of 20 smaller than the 
ls2s^2p(^P) ls^2s'^{}S) rate, so it only increases the fluorescence yield by a few percent 
at low Z. As Z increases, however, eventually even this reduced radiative rate dominates 
the autoionization rate, giving 

1 . 



(9) 



(10) 



-6- 



The spin-orbit interaction also affects tfie computed fluorescence yield, primarily by 
mixmer the ^Pi and ^Pi levels. The mixing fraction, while only about 6.3x10 ^ at Z = 5, 
has a dependence, and eventually becomes quite significant, reaching 0.117 at Z = 30. 
As a result, the "^P-j^" level (this is now just a label used to indicate the dominant term of 
a level) has an increased fluorescence yield, and we get that the intermediate coupling (IC), 
configuration-averaged fluorescence yield, including CI, behaves as 

> Clsci. (11) 

Thus we see that CI and the spin-orbit interaction each cause an increase in the computed 
fluorescence yield as Z is increased. 



2.2. Earlier Be-like Fluorescence Data 

The approach of Kaastra & Mcwe (1993) for this particular Be-like series was to ne- 
glect spin-orbit and CI effects, and to assume that the hydrogenic approximation is valid 
throughout the series. Furthermore, they used configuration-averaged values for the B II au- 
toionization rate and emission oscillator strength, which were computed by McGuire (1969), 
and the experimental values of u from Lotz (1967, 1968), to obtain the ratio Ar/{Ar + Aa) 
required for determining ^ using Eq. 7. This is not the same as the desired configuration- 
averaged fluorescence yield ^lssc in Eq. 9 - the ratio of the averages does not equal the 
average of the ratios: 



A,(2'5+lp) 



A,(25+1P) + A„(25+1P) 



12) 



We first address the accuracy of the computed autoionization rates and emission os- 
cillator strengths in the next subsection, and then address the validity of the hydrogenic 
approximation in the following subsection. Fluorescence yields are presented in the last sub- 
section, where the incorrect averaging and neglect of CI and spin-orbit effects by Kaastra & 
Mewe (1993) are addressed. 



2.3. Atomic Calculations for B II 



In order to calculate the transition matrix elements appearing in the expressions for 
the radiative and autoionization rates (Cowan 1981), it is first necessary to produce atomic 
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wave functions. McGuire (1969) used the Herman-Skillman approximation in determining 
the (single-configuration) wave functions, whereby all electrons (i.e., the Is, 2s, 2p, and 
continuum ones) arc cigenf unctions of a common central potential; as stated by McGuire 
(1969), this "neglect (s) ... exchange and correlation effects." Furthermore, this potential 
rV{r) is approximated by "a series of straight lines" in order to yield piece-by-piece analytic 
Whitakker functions. Here we are concerned with the validity of these approximations, 
given that more rigorous calculations can be easily performed using today's state-of-the-art 
technologies. 

For the present study, we use the program AUTOSTRUCTURE (Badncll 1986), which 
generates Slater-type Is, 2s, 2p, and distorted-wave continuum orbitals. In order to compare 
with the results of McGuire (1969) for B II, and with Kaastra & Mewe (1993) as we scale 
from Z — 5 to Z — 30, we first performed single-configuration LS calculations. For the more 
rigorous calculations that we compare to other theoretical results and that we recommend as 
the definitive data, we also included CI - ls2s^2p -|- ls2p^ for the inner-shell vacancy state 
and ls^2s^ + ls^2p^ for the final radiative decay state - and spin-orbit effects. The two 
accessible continua were described as ls^2sep and ls^2pes, where el denotes a continuum 
distorted wave. 

Given atomic wave functions, McGuire (1969) computed the configuration average (CA) 
radiative and partial autoionization rates in Eqs. 1-3. The emission oscillator strength given 
is thus 



whereas for the total autoionization rate, the CA rates for the processes in Eqs. 2-3 were 
used, that is. 



f{CA) = + ^/fp) 



= 0.0377 , 



(13) 



Aa{CA) 



= Aal{CA)+Aa2{CA) 

= 2.37 X 10"^ a.u. , 



(14) 



where 



Aal{CA) 



\a^,{'p) + \a^,{'p) 




+-{27r[i?o(ls,ep,2s,2p)f} . 



(15) 
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and 

Aa2{.CA) = -^Aa2{^P)+'^-Aa2{^P) 

= 27r[i?o(ls,es,2s,2s)]^ , (16) 

since our calculations indicate that Aa2{^P) — A^^i^P)- Here R\{nili,n2l2,n3ls,nj4) is a 
Slater integral of multipole A (Cowan 1981), and el represents the outgoing Z-wave continuum 
electron orbital. (The expressions in Eqs. 15 and 16 are equivalent to those in Eq. 6 
of McGuirc (1969) for inequivalent electrons and single-p orbital occupation, considering 
the different continuum normalization used by McGuire (1967)). Note that the partial rate 
AaiC'P) in Eq. 15 is greatly suppressed relative to the Aai{^P) rate due to a near cancellation 
of monopole and dipole Slater integrals. (Indeed, it was due to this near cancellation of Slater 
integrals that Caldwell et al. (1990) explained why the inner-shell photoexcited ls2s^2p(^P) 
resonance in Be I preferentially decayed - by two orders of magnitude - to the ls^2p(^P) +e~ 
channel, compared to the ls^2s(^5') + e~ channel.) Thus the configuration average partial 
rate Aai{CA) will be larger than the partial rate Aai{^P), and hence the configuration 
average total rate Aa{CA) will be larger than Aa{^P). 

Since we are interested in computing ^(^-P), which requires Aa(^P) and f^P), we have 
converted the reported values from McGuire (1969) to the ^P values (the Slater integrals 
were also given in that work) . We get the following values 

/(^P, McGuire) = 0.1508 (17) 
A„(^P, McGuire) = 1.692 x 10"^ a.u. , (18) 

which compare fairly well with our results obtained using AUTOSTRUCTURE: 

/(^P, present) = 0.1519 (19) 
A„(^P, present) = 1.045 x 10"^ a.u. (20) 



In summary, we find that the earlier results for B II of McGuire (1969) are consistent 
with ours. However, for the astrophysical plasma modeling purposes we have in mind, one 
really requires the configuration average fluorescence yield, not the ratio of the averaged 
radiative and total rates used by Kaastra & Mewe (1993) 

AriCA) 



Ar{CA)+A,{CA) 

1 , (21) 



due to A^ and Aa scaling as and respectively, in the hydrogenic approximation. Equa- 
tion 21 differs from the correct ^lssc given in Eq. 9. First, we have Aa{CA) > Aa{^P) due to 
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the near cancellation in the 2s2p — > Isep partial autoionization rate, so at low Z, where 
Ar <^ Aa, we have ^{K^M) < ^lssc- Second, when CI and spin-orbit effects are ignored, 
as they were in Kaastra & Mewe (1993), the fluorescence yields differ asymptotically by a 
factor of 4, 

lim i^- L = 4 , (22) 

1;lsSC 

as can be seen by comparing Eqs. 9 and 21. Of course, CI needs to be included for all Z, 
whereas spin-orbit mixing needs to be included at higher Z, and both ^{KSzM) — > 1 and 
Cicci — > 1 as Z — > oo. However, in the intermediate Z range, it can be shown that the 
Kaastra & Mewe (1993) results are still larger than the ICCI results. 



2.4. Validity of the Hydrogenic Approximation 

In order to assess the vahdity of scahng the B II results along the isoelectronic series, 
we computed both the autoionization rate Aa and emission oscillator strength / for all 
Be-hke ions up through zinc, first neglecting spin-orbit effects. In Fig. 1, it is seen that 
neither of the two is independent of the nuclear charge Z at the lowest stages of ionization 
- the emission oscillator strength increases by about 2/3 in going toward the highly- ionized 
regime whereas the autoionization rate more than doubles. Furthermore, by choosing the 
scale so that our two quantities coincide for B II, it is seen that the important ratio Aa/f 
appearing in Eq. 7 increases by roughly 25% by the time Zn XXVII is reached. Thus 
the assumption of pure hydrogenic scaling by Kaastra & Mewe (1993) alone introduces an 
uncertainty at the highly-charged end of this series. Due to the stronger Z dependence at 
the near-neutral end, together with the greater sensitivity to the atomic basis used in this 
region, we recommend that if scaling along an isoelectronic sequence is to be performed, the 
better starting point would be at the highest Z desired, extrapolating the rates to lower Z 
members. Of course, given the ease of determining atomic rates with modern computing 
capabilities, the most reliable approach is to calculate the fluorescence yield directly rather 
than resort to questionable scaling methods. 



2.5. Fluorescence Yield Results 

While the assumption of hydrogenic scaling introduces an ^ 25% inaccuracy in Aa/f, 
the initial quantity being scaled in Kaastra & Mewe (1993) - the ratio of averages rather 
than the average of ratios - is really not the desired quantity to be scaled in the first place. 
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Together, these approximations lead to an uncertain prediction for the fluorescence yield. 
In Fig. 2 (and Table 1), we compare various results for ^ along the Be-like sequence, where 
it can be seen that our single-configuration LS results differ greatly from those of Kaastra 
& Mewe (1993), especially at higher Z] here, especially, their results are expected to differ 
from the correct single-configuration values due to their incorrect asymptotic value given 
by Eq. 21. A more disturbing result was found when we tried to repeat their calculations, 
i.e., when we used Eq. 7, with the ratio of Aa{CA)lf{CA) taken from McGuire (1969), and 
the energies lo taken from Lotz (1967, 1968). Whereas these scaled results exhibit a smooth 
monotonic increase with nuclear charge Z, those of Kaastra & Mewe (1993) are somewhat 
irregular, showing unphysical dips, and do not agree with what we tried to reproduce, given 
their stated method. Either way, the results of Kaastra & Mewe (1993), or our scaled ones 
using the B II results of McGuire (1969), initially underestimate our results at lower Z, and 
then overestimate our (LSSC) results by almost a factor of 3 for the highest Z = 30. 

To our knowledge, there have been two other calculations for the fluorescence yields of 
some members of the Be-like sequence: those of Behar & Netzer (2002) using the HULLAC 
codes (Bar-Shalom et al. 2001) and those of Chen (1985) using a multiconfiguration Dirac- 
Fock (MCDF) method. In both cases, CI and spin-orbit effects were included. Here we do 
the same, first adding the important 2s^ — > 2p^ CI discussed earlier to the LS calculations in 
order to see that this effect increases the Z = 30 fluorescence yield by about 30%. Then when 
spin-orbit effects (and other higher-order, relativistic effects) are included in our intermediate 
couphng calculation, there is a further increase in the fluorescence yield by about 20% 
more. In comparison with the other two calculations along this series, there is overall good 
agreement with these IC results. 



We turn now to the simplest closed- (outer) shell case of a Is- vacancy in F-hke ions, 
giving the ls2s^2p^(^S) state which decays as 



3. Case Study of the F-Like Fluorescence Yields 




( W2p\^S)es 
I W2s2p^{^^^P)ep 
\ ls^2s^2p\^P, ^D, ^S)es,ed 



(23) 




(24) 



Again, only one photon, or one electron, can be emitted, which simplifies the analysis consid- 
erably (when spin-orbit effects are considered, the final ionic term in Eq. 23 is fine structure 
split into the ground ls^2s^2p^(^P3/2) level and the metastable ls^2s^2p^(^Pi/2) level). Since 
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this is a closed-shell system, the Herman- Skillman method for the important 2p electrons 
is expected to be more accurate than for B II. Indeed, as stated by McGuire (1969), "in 
stripping away electrons (in reducing to a closed-shell system), ... we should be increasing 
the applicability of the common central-field approximation." Furthermore, there is only 
one ls2s^2p^ Is-vacancy state, rather than the two wc had for the Be-like sequence, and no 
other intrashell configurations to CI mix with, so we do not need to consider population of 
non-fluorescing states by CI or spin-orbit mixing, nor do we have to consider configuration 
averaging issues. Consequently, a single configuration LS coupling calculation is sufficient to 
determine accurate Aa, /, and ^ values for the term. 

As a result, the computed values of the autoionization rate and emission oscillator 
strength given by McGuire (1969) agree quite well with our values, as seen in Fig. 3 and 
Table 2. However, both of these values depend on the internuclear charge Z, giving a ratio 
Aa/ f that increases by about a factor of 1/2 in going from Ne II to Zn XXII. Thus the 
scaled fluorescence yield ^, using Eq. 7, the ratio Aa/ f from McGuire (1969), and u from 
Lotz (1967, 1968), increases relative to the actual computed value, as is seen in Fig. 4 and 
Table 2. The more troublesome news in this flgure is the actual tabulated values of Kaastra 
& Mewc (1993) - their values do not follow our attempt at reproducing those results, but 
rather tend to follow our computed values, except for certain unphysical dips. Nevertheless, 
the results reported by Kaastra & Mewe (1993) for F-like ions are not plagued by as many 
uncertainties as those for Be-hke ions. We also see in Fig. 4 that the HULLAC results are 
in good agreement with our present ones (the results reported earlier by Behar & Netzer 
(2002) only considered fluorescence into the ls^2s^2p^(^P3/2) level, which includes only 4 of 
all 6 magnetic sublevels of the ls^2s^2p^(^P) conflguration; therefore, those values must be 
multiplied by about 3/2 to account for fluorescence into the two ls^2s^2p^(^Pi/2) sublevels 
as well. Furthermore, the earlier HULLAC result for F+ was erroneously listed incorrectly, 
and here we have given the actual computed value that should have appeared). 



4. Summary and Conclusion 

The inaccuracies we have discovered in the reported results of Kaastra & Mewe (1993) 
for Be-like ions are as follows: 

1. The computed atomic data for B II are used in the form Aa{C A) / f (C A) , that is, the 
radiative and autoionization rates have been averaged over the and ^P configura- 
tions, whereas the desired quantity for plasma modeling applications is ^icci and is not 
the same thing, differing qualitatively and quantitatively, especially in the asymptotic 
high-Z limit. 
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2. The hydrogenic scaling assumed is invalid. The autoionization rates, the emission 
oscillator strengths, and even the ratio A^j j are not independent of nuclear charge Z. 

3. The tabulated data of Kaastra & Mewe (1993) do not seem to follow the results we 
obtain when we try to reproduce their stated method using Eq. 7, with AJyCA) / f{CA) 
from McGuire (1969), and uj from Lotz (1967, 1968). 

4. The calculations of Kaastra & Mewe (1993) neglected CI and spin-orbit effects as they 
scaled to higher Z. 

For F-hke ions, items 1 and 4 are not issues since there is only one inner-shell vacancy 
term. However, points 2 and 3 still apply for the F-like sequence. For plasma modeling 
purposes, we recommend our ^icci for the Be-like sequence and our ^ for the F-like sequence. 

In conclusion, we propose that, given the many uncertainties discovered, the entire data 
base of Kaastra & Mewe (1993) should be reevaluated. While we have focused on systems 
that can emit only one photon or one electron, their comprehensive tabulation also includes 
data for ions with n > 3 shells occupied; these can emit multiple electrons and/or photons 
through numerous cascading channels, compounding the inaccuracies we have discovered. 
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Table 1. Emission oscillator strengths, autoionization rates, fluorescence yields, and 

photon energies for Be-like ls2s^2p ions. 
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0.1519 
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0.1045 
0.1692'" 


0.0014 


0.0006 


6.751 


0.0006 




0.0011 


0.0011 
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0.1712 


0.1194 


0.0032 


0.0019 


10.349 


0.0013 


0.0024 


0.0025 


0.0025 
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0.1859 


0.1328 


0.0061 


0.0052 


14.721 


0.0027 


0.0045 


0.0048 


0.0048 


8 


0.1972 


0.1442 


0.0106 


0.0096 


19.866 


0.0049 


0.0079 


0.0083 


0.0083 
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0.2062 


0.1540 


0.0168 


0.0154 


25.753 


0.0081 
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0.1369 


0.1671 


112.664 


0.1357 


0.1273 0.1237 
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0.3838 
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0.3810 
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0.4200 


0.2394 0.2633 
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27 
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0.2169 
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0.5267 


261.659 


0.4584 




0.2585 


0.2916 


28 
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0.2182 
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0.5836 


282.129 


0.4960 




0.2712 
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29 
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0.6215 


303.333 


0.5322 




0.2840 


0.3339 


30 


0.2566 


0.2205 


0.2261 


0.6322 


325.310 


0.5668 




0.2969 


0.3553 



^Present LS results for emission from the term (dimensionless). 

Present LS results, autoionization from the term (in units of 10~ a.u., 1 a.u.= 4.13 X 10^^ s-l). 
"^Present LS results using a single configuration, one fourth the term fluorescence yield (dimensionless). 
''Kaastra & Mewe (1993). 

«Lotz (1967, 1968) (in a.u., 1 a.u.= 27.211 eV). 

'Obtained using Eq. 7 with Aa(CA)//(CA) for B II from McGuire (1969) and ixi from Lotz (1967, 1968). 
SBehar & Netzer (2002), averaged over the ^Pi and ^Po,i,2 levels. 
''Chen (1985), averaged over the ^P\ and ^Po,i,2 levels. 

'Present LS results, including configuration interaction (CI), averaged over the ^P and '*P terms. 
jPrescnt intermediate-coupling (IC) results, including CI, averaged over the ^P\ and ^Po,i,2 levels. 
■"McGuire (1969). 
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Table 2. Emission oscillator strengths, autoionization rates, fluorescence yields, and 
photon energies for F-hke ls2s^2p^(^S) ions. 
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0.3429 
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0.2751 
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0.4041 
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0.4531 
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284.040 
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0.4734 


305.465 


0.6230 




30 


0.3406 


0.2420 


0.5207 


0.4758 
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0.6554 





'^Present results (dimensionless). 

Present results (in units of 10~ a.u., one a.u.= 4.13 X lO^s s-i). 
"^Kaastra & Mewe (1993). 

'^Lotz (1967, 1968) (in a.u., 1 a.u.= 27.211 eV). 

"Obtained using Eq. 7 witli Aa(,CA)/f{CA) for B II from McGuire (1969) 

and UJ from Lotz (1967, 1968). 

'Behar & Netzer (2002). 
SMcGuire (1969). 
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^ A^(McGuire, units of 6.88 x 10" a.u.) 

(McGuire, dimensionless) 
^ Present (units of 6.88 x 10"^ a.u.) 
o Present f (dimensionless) 
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Fig. 1. — Present LS autoionization rates Aa (in units of 6.88 x 10~^ a.u., open squares) and 
emission oscillator strengths / (dimensionless, open circles) for Be-like ls2s'^2p(^P) ions as 
a function of the nuclear charge Z. The autoionization rate and emission oscillator strength 
from McGuire (1969) for B II are shown by the solid square and circle, respectively. 
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Fig. 2. — (a) Comparison of various computed and inferred fluorescence yields ^ for Be- 
like ls2s^2p ions: present LS results in the single-configuration (SC) approximation - solid 
diamonds; present LS results with configuration interaction (CI) included - solid squares; 
present intermediate couphng (IC) results with CI included - sohd circles; HULLAC results 
from Behar & Netzer (2002) - open squares; MCDF results from Chen (1985) - open dia- 
monds; Kaastra & Mewe (1993) - open circles; results when we scale McGuire (1969) B II 
results, using Eq. 7 and the uj from Lotz (1967, 1968) - open triangles, (b) Same as (a) 
focusing on the \ow-Z region. 
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□ □ Present (units of 4.89 x 10" a.u.) 

o o Present f (dimensionless) 
• f (McGuire, dimensionless) 
■ A„ (IVIcGuire, units of 4.89 x 10"^ a.u.) 
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Fig. 3. — Present LS emission oscillator strengths / (dimensionless, open circles) and au- 
toionization rates Aa (in units of 4.89 x 10^^ a.u., open squares) for F-like ls2s'^2p^(^S) ions 
as a function of the nuclear charge Z. The emission oscillator strength and autoionization 
rate from McGuire (1969) for Ne II, as used by Kaastra & Mewe (1993), are shown by the 
solid circle and square, respectively. 
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Fig. 4. — (a) Comparison of various computed and inferred fluorescence yields ^ for F-like 
ls2s'^2p^{^S) ions: present LS results - solid circles; Kaastra & Mewe (1993) - open circles; 
results when we scale McGuire (1969) Ne II results using Eq. 7 and the uj from Lotz (1967, 
1968) - open triangles; HULLAC results of Behar & Netzer (2002) - open squares, (b) Same 
as (a) focusing on the \ow-Z region. 
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